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Abstract. In this paper we study finite dynamical systems with n functions 
acting on the same set X, and probabiUties assigned to these functions, that it 
is called Probabilistic Regulatory Gene Networks (PRN) in S; . This concept 
is the same or a natural generalization of the concept Probabilistic Boolean 
Networks (PBN), introduced by I. Shmulevich, E. Dougherty, and W. Zhang 
in particularly the model PBN has been using to describe genetic networks 
and has therapeutic applications, see 1101 . In PRNs the most important ques- 
tion is to describe the steady states of the systems, so in this paper we pay 
attention to the idea of transforming a network to another without lost all the 
properties, in particular the probability distribution. Following this objective 
we develop the concepts of homomorphism and e-homomorphism of proba- 
bilistic regulatory networks, since these concepts bring the properties from 
one networks to another. Projections are special homomorphisms, and hey al- 
ways induce invariant subnetworks that contain all cycles and steady states in 
the network . 



Introduction 

Genes can be understanding in their complexity behavior using models according 
with their discrete or continuous action. Developing computational tools permits 
describe gene functions and understand the mechanism of regulation 0] |nj . This 
understanding will have a significant impact on the development of techniques for 
drugs testing and therapeutic intervention for treating human diseases [51 151 [TU] . 

We focus our attention in the discrete structure of genetic regulatory networks, 
instead of, its dual moving continuo-discrete. Probabilistic Gene Regulatory Net- 
work(PRgN) is a natural generalizations of the model Probabilistic Boolean Net- 
work (PBN), introduced by I. Shmulevich, E. Dougherty, and W. Zhang in 'W . The 
mathematical background of the model PgRN, is introduced here, for simplicity 
we work with functions defined over a set X to itself, with probabilities assigned 
to these functions. X is a set of states of genes, for example X = {0, 1}", if our 
network is a Boolean network. Working in this way, we can observe the dynamic 
of the network indeed focus our attention in the description of functions. The set 
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X can be a subset of {0, 1}", and we can extend some classical ideas to regulatory 
network, such as invariant subnetworks, automorphisms group, etc. In particular 
if X is a vector space over a finite field, the functions are lineal functions, then 
we can use linear Algebra to describe the state space. Mapping are important 
in the study of networks, because they permit to recognize subnetworks, in par- 
ticular determine when two networks are similar or equivalent. Special mappings 
are homomorphisms and e-homomorphisms, we use both to describe subnetworks 
and similar networks. An homomorphism transform a network to another in such 
a way the discrete structure giving by the first network can lives in part of the 
other one, or these two networks are very similar but no equals, in particular in 
the probabilistic way. An e-homomorphism is the same but with the condition that 
the probability distributions of the networks are close, and we use a preestablishes 
< e < 1 as a distance between the probabilities. 

1. Finite dynamical systems and probabilistic Boolean networks 

Two finite dynamical systems (X, /) and {Y, g) are isomorphic (or equivalents) 
if there exists a bijection (f) : X such that (j> o f = g o cj)^ (or/ = (/)~^ ° g ° 4>)- 
If (j) is not a bijection map then cj) is an homomorphism. 

If y C X is such that f{Y) C Y then (y, f\y) is a sub-FDS of [X, /), where /|y 
is the map restricted to Y. There exists naturally an injective morphism from Y to 
X called inclusion and denoted by l. The state space of a FDS {X, /), is a digraph 
whit vertices the set X , and with an arrow from u to v ii f{u) — v. 

Example of FDS-homomorphism The FDSs X = ({0, 1}^, y) = {xy,y)), 
and Y = ({0, 1}^, /2(a;, y) = {x,{x + l)y)) are isomorphics, because their state 
spaces are isomorphics. 

(1.0) -^f^ (0,0) (0,1) 

(1,1) 

(1.1) (1,0) (0,0) 

(0,1) 

In fact, the isomorphism (p : {0,1}^ {0,1}^ is the bijection 0(1,0) = (1,1), 
0(0,0) = (1,0), 0(0,1) = (0,0), and 0(1,1) = (0,1). The following is an example 
of homomorphism (inclusion) with Z = {{(0, 0), (1, 0)}, /i}. 

(1,0) (1,0) (0,1) 

(0,0) O (0,0) O (1,1) 

A Probabilistic Boolean Network A ~ (V, F, C) is defined by the following sort 
(type) of objects [5| ?]: a set of nodes (genes) V = {xi, . . . ,a;„}, Xi e {0, 1}, for 
all I- a family F = {i^i, F2, . . . , F„} of ordered sets = {/f \ , . . . , of 

Boolean functions /j*^ : {0, 1}" {0, 1}, for all j called predictors; and a Hst 
C ~ (Ci, . . . , C„), Ci = {ci \ . . . , c^(^j)}, of selection probabilities. The selection 
probability that the function /j'' is used for the vertex i is c^*' — Pr{f^'^^ = fj^^}- 
The dynamic of the PBN is given by a vector of functions ffe = ifll\fj^\ • ■ • , /fc"^) 
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for 1 < ki < and fl 6 F,, where k = [ki, . . ., kn], 1 < ki < The map : 
{0, 1}" {0, 1}" acts as a transition function. Each variable Xi G {0, 1}" represents 
the state of the vertex i. All functions arc updated synchronously. At every time 
step, one of the functions is selected randomly from the set Fi according to a 
predefined probability distribution. The selection probability that the transition 
function fj. = {.fi]^^,.fi^J: ■ ■ ■ , fk^) is used to go from the state u G {0, 1} to another 
state ffe(u) = G {0, 1}" is given by 

n 
i=l 

The dynamical transition structure of a PBN can be described by a Markov chain 
with fixed transition probabilities. Tlicrc^ arc two digraphs structures associated 
with a PBN: the low-level digraph F, consisting of genes functions essentiality 
relations; and the high-level digraph which consists of the states of the system and 
the transitions between states. The matrix T associated to the high level digraph 
formed by placing p(u,v) in row u and column v, where u, v G {0, 1}" is called 
the transition probability matrix or chain matrix, p{u, v) = J2t^\t^{u)=v "-ffc • 

2. Probabilistic Regulatory Gene Networks 

A Probabilistic Gene Regulatory Network (PRN) is a triple X = {X, F, C) where 
X is a finite set and F = {/i, . . . , /„} is a set of functions from X into itself, with 
a list C = (ci, . . . ,c„) of selection probabilities, where Cj = p{fi)- We associate 
with each PRN a weighted digraph, whose vertices are the elements of X, and 
if u,w G X, there is an arrow going from u to w for each function fi such that 
fi{u) = V, and the probability c,; is assigned to this arrow. This weighted digraph 
will be called the state space of X. In this paper, we use the notation PRN for one 
or more networks. 



Example 2.1. 

If X = {0, 1}2, F = y) = {x, y), f2{x, y) = (x, 0), 

h{x,y) = {l,y),f4{x,y) = (1,0)}; and C = {.46, .21, .22, .11}, the state space of 
X = {X, F, C) is the following: 



.67 

O (0,0) 

•33} 

0(1,0)^(1,1)0 



21 (0,1) 0-46 

^.22 J. ^ 

68 



.11 J 



.67 





.33 





.21 


.46 


.11 


.22 








1 











.32 


.68 



3. HOMOMORPHISMS AND e-HOMOMORPHISMS OF PRN 

If C is a set of selection probabilities we denote by x the characteristic function 
over C. That is x : {0} ^ {0, 1} such that x(c) = 1, if c ^ and x(0) = 0. Let 
X^ = {Xi,F= (/i)r=i, C) and X2 = {X2,G= igj)JLi,D) be two PRN. 



Definition 3.1 ( Homomorphisms of PRN). A map cf) : X\ ^ X2 is an homo- 
morphism from X\ to X2, if for all fi there exists a gj, such that for all u, v in 

(1) (j) o fi = Qj o (j); and (2) xidg^ (^(w), (l>{v))) > x{cfi {u, v)). 
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X2 — ^ X2 

If (j) : Xi ^ X2 is a bijective map, then ^ is an isomorphism . 

Example 3.2 (PRN-Homomorphism). 

lfX= {X; F; C) is the PRN in Example 1, and Xi = {X; F' = {/i, /2, /a}; C" = 
{.47, .28, .25}) is a new PRN over the same set X with different probabihties and 
only three functions. 



X 



47 



0(0,0) ^-28 (0,1) "0 
Qi (1,0)^(1,1)0-72 



.46 

21 (0,1) O 



0(0,0) 

i /.lli-22 

01 (1,0)^(1,1)0- 



.75 .25 

.28 .47 .25 

10 

.28 .72 



67 





.33 





21 


.46 


.11 


.22 








1 











.32 


.68 



The homomorphism (j) : Xi X is a. bijective map, (^(a;) = x, over the set of states, 
but an inclusion over the set of arrows, because the arrow going from (0, 1) to (1,1) 
in X doesn't appear in Xi. The first condition for homomorphism is obvious. The 
condition (2) holds, because the inclusion of arrows. The two transition matrices 
are connected by this inclusion, since if the place ij in the first matrix ^ then this 
place is 7^ in the second network too. The two PRN are not isomorphics because 
the probabilities are not equals. Since, there are no specific condition about the 
probability distribution in both PRN, we include a third condition, obtaining in 
this way a new concept that we will call e-homomorphism of PRN. 

Condition (3) for e-Homomorphism The distributions of probabilities fol- 
lowing the homomorphism are enough close. An e- homomorphism is an homomor- 
phism that satisfies the condition, for all i, j, max\p{ui,Uj) — p{(j){ui) , (l>{uj))\ < e, 
where e > is a real number that we previously determine for the applications. 

As a consequence of this condition, if we use a test, as Kolmogorov-Smirnov test, 
the difi^crences between the two distributions are < e again. In order to determine 
e for the homomorphism, we use the transition matrices. In the above example 
e= .11. 



Ti-T = 



.08 





-.08 





.07 


.01 


-.11 


.03 




















-.04 


.04 



Conclusion 

If the homomorphism is a bijective map like here, the transition matrices Ti and 
T2 have the same order, and X^"=i(T'i — T2)ij = 0, for j = l,n 



Theorem 3.3. If(p:Xi^X2 is an e-homomorphism, then the transition matrices 
Ti and satisfy the condition: 
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max\{Ti^'')ij — (Tfi/^)ij\ < e, for all possible i and j, and all n > 1. // the homomor- 
phism is injective and e < 1, the steady state of Ti and the steady state ofTfj, are 
close, that is satisfy 
ki - TT^I = maXi\Tri{i) - TT^{i)\ < e . 

Proof. It is clear if we do the following 

\p{uj\u))~p{cl,{u),cl,{f^{u))\ = 
\p{u, f{u))p{f{u), f\u)) - p{c^{u), cj,{f{u)))p{4>{f{u)),^{f\u)))\ 
< \pifiu),fiu))\\piu, fiu)) - Hfiu)))\ + 

\p{q^iu),c^ifiu)mpifiu),fiu))~pmiu)),^if\u)))\ < e 
Then our aim holds. □ 

4. Algebra of Probabilistic Regulatory Networks 

Sum of two PRN 

Let Xi = {Xi,F = C) and X2 = (X2, G = {gj)JLi,D) be two PRN. The 

sum Xi®X2 = {X1UX2, F V G, C V D) is a PRN where 

(1) Xi(jX2 is the disjoint union of Xi and X2. 

(2) the function hij = [fi V gj) is defined by hij{x) — fi{x) if a; G Xi and 
hij{x) = gj{x) ii X ^ X2- 

(3) the probability ) — aVdj, that is p{hij) — a if hij — fi or pihij) = dj 
if hij = gj. 

If Ti and T2 are the transition matrices of Xi and X2 respectively, Then T = 
T \ 

„ I is the transition matrix of Xi (B X2. 
U I2 J 

Example 4.1. 

An example of sum is the PRN obtained by summing the same PRN twice, 
X ® X. To make the disjoint union, we subindicate X with for the first X and 
with 1 for the second X. That is, the new set is 

XoUXi = {(0, 0, 0), (0, 1, 0), (0, 1, 0), (1, 1, 0)} 

U{(0,0,1), (0,1,1), (0,1,1), (1,1,1)}. 

The digraph is: 

(1.1.0) (1,0,0) (0,0,0) 

(0,1,0) 

(1.1.1) (1,0,1) (0,0,1) 

(0,1,1) 

This is a way to construct a PRN over {0, 1}" using either one or two PRN over 
{0, l}"-\ since 2"-^ + 2"-^ = 2". 

Superposition 

It is clear that a PRN is the superposition of several Finite dynamical Systems 
(FDS)in] over the same set X with probabilities assigned to each FDS. Since each 
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functions defined over a finite field can be wrote as a polynomial function, we will 
use this notation for functions over a finite field, 0. If X = {0, f } = Z2, the finite 
field of two elements, all the FDSs over X have one of the following state space, 
where fi{x) = x; f2{x) = 1; /3(a;) = 0; fi{x) = x + l,\/x & X: 





L2 


L3 


Li 
















i 


T 


n 


1 


1 


1 


1 



If Pi denotes de probability assigned to Li, and Ti denotes its transition matrix, 
then the set of all PRN is described as follows. 

^ \ P3+P4 P1+P2 J ^ 

We denote by L1L2 the superposition of Li and L2, and similarly L1L3 is the 
superposition of Li and L3. The state spaces are the following: 



L1L2 


L1L3 


L1L4 










IP' ' 




P4 

n 


1 0^ 


I 0^*^ 


I OP' 



L2L3 -/j2-/j4 

0^*^ 

|P2-[-P3 |P4j-P2+P4 ' 

1 O^'^ 1 0^^ 

L3L4 

OP" 
p* , > 

j jP4+P3 
1 

For example, with transition matrices 

r. = T. + r. = (« «)r„ = T. + T,= (;^ ;J 

Product of tw^o PRN 

Let Xi = {Xi,F= ifiYl^^.C) and X2 = [X2, G = {gj)JLi,D) be two PRN. The 
product A:"! X A:'2 = (Xi X X2, F X G, C A D) is a PRN where 

(1) Xi X X2 is the cartesian product of Xi and X2- 

(2) the function hij — {fi,gj) is defined by 

hij{xi,X2) = {fi{xi),gj{x2)) 

for xi G Xi, and X2 G X2. 

(3) the probability is a function of Cj and dj, for example p{hij) — '^'"^''^ . 

Example 4.2. 
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The product L1L2 x L1L3 is the PRN with four states {(0,0), (0, 1), (1,0), (1, 1)} 
and four functions 

fii{x,y) = {x,y),fi3(x,y) = (a;,0), 
f2i{x,y) = (l,y), f23{x,y) = (1,0). 
The state space is the following: 

L1L2 X L1L3 

P11+P13 Q (0^0) ^Pi^ (0,1) O^'" 

P23+P21 I / IP21 
-I- 1/ P23 -1- 

(1,0) < — (1,1) op^'+p^' 

P13+P23 

The transition matrix is the following 

/ Pii + Pi3 P23 +P21 \ 

rp _ P13 Pu P23 P21 

1 

\ P13+P23 P11+P21 J 

4.1. Linear Probabilistic Regulatory Networks. A linear PRN is a superpo- 
sition of linear FDS. A linear FDS is a pair {X, /) where / is a linear function, and 

X is a vector space over a finite field. So, a linear PRN is a triple {X, (/j)™ 1, C), 
where X is a finite vector space, the functions fi'.X^X are linear functions, and 
C = {ci = p{fi)}- The set X has cardinality a power of a prime number and each 
linear function is determined by its characteristic polynomial and the companion 
matrix. 

If X = Z3 = {0, 1, 2} is the field of integer modulo 3, then the linear functions 
are: fi{x) = x, f2{x) = 2x, and fsix) = for all x G Z3. So, the linear PRN are 
the following: 

{/l,/2} {/l,/3} 



(jPi 1 ^P2 2 1 Op, Op, 2 

{/2,/3} {/l,/2,/3} 



yp3'\^ yp3'\^ 

1 P2 2 Opi 1 2 Opi 
If X = Z2 X Z2 is the vector space with 4 elements over the field Z2, then there 

are 4 linear FDS not isomorphics. In fact, using matrix, the possible characteristics 
polynomials p/(A) are: A^, A-^ + A + 1, A^ + A + 1. The companion matrices of 
these linear functions are: 

-.(0 s)-Mo DHi ;)'^'(? 

Then the FDS associated to this matrices are: 



Ai A2 

0(0,0) ^ (1,0) 0(0,0) ^ (1,0) 

T \ 

(0,1) (1,1) 0(0,1) ^ (1,1) 
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P3 



O(0,0)fi(l,0)O 
Pi I\P1 Pi -I-^PI 



A3 A4 

0(0,0) 0(1,0) 0(0,0) (1,0) 

/ T 

0(0,1) 0(1,1) (0,1) ^ (1,1) 
The linear PRN with two functions are the following: 

Ai,A2 AuAs 

0(0,0)^(1,0) oi 

pi I'^pi 

0(0,1)^(1,1) - O (0,1) (1,1) of^ 

0^0,0)^(1,0) 0^ (0,0)<^PMl,0)Of=* 

Pl I / |P4 

(o,i)fi(i,i) oMo,i) (1,1)0- 
oMo,o) ^''^ (1,0) 0^0,0) (1,0)0- 

P4 ^I^P* P<1 ^fP-* 

O- (0, 1) (1, 1) - O (0, 1) -^P^ (1, 1) o^*^ 



'3 



5. Invariant Subnetworks and Projections 

A subnetwork Y Q X oi X = {X, F, C) is an invariant subnetwork or a sub- 
PRN of X if fi{u) € Y for all u&Y, and fi S F. Sub-PRNs are sections of a PRN, 
where there aren't arrows going out. The complete network X, and any cyclic state 

with probability 1, arc sub-PRNs. An invariant subnetwork is irreducible if doesn't 
have a proper invariant subnetwork. An endomorphism is a projection if tt^ = tt. 



Theorem 5.1. If there exists a projection from X to a subnetwork y then y is an 
invariant subnetwork of X. 

Proof. Suppose that there exists a projection tt : X ^ Y. If y G F, by definition 
of projection 7r(j/) = y, and /i(7r(j/)) = n{gj{y)). Therefore all arrows in the 
subnetwork Y are going inside Y, and the network is invariant. □ 

Example 5.2. 

The PRN X has two invariant subnetworks with projections 
7ri(a;,y) = (a;,0) and7r2(a;,y) = (l,y). 

.67 .46 

O (0,0) ^-21 (0,1) O 

Q9 .68 
Ql (1,0)^(1,1) O 



.67 .67 

00 O (0,0) ^ ^1 X 

1 1 



-Q >J \\J,\J) „ 

'^l .33 I ~ .33 I "^l ' TTs 



^1 ^(I'O) 5,0^1,0)^(1,1)0 

— 1 32 -68 

52 O^ ^ 1 
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Checking the probabihties for tti and TT2, we have ei — .68; and 62 
We can observe that X = Si x 82- 



.67. 



Example 5.3. 

The subnetwork Xi 
{{0,iy, F,C). 



{{{x,y,l)}, F,C) is an invariant subnetwork of X 



000 



.549 

^451 



X 



100 

.995 
i 

010 

.378 

o 



.005 



.622 



001 



110 



Xi 

.113, .456 

.544 

\ 
.002 



,998 



.448 

101 o 

.439 

i 

.337 

Oil c 

.663.011 

IT 

.989 

111 c 



00 



.1 13, .45 6 
.544 



X-[ = Xi 



10 

.439 

i 

01 

.011.663 

Ti 
11 



.448 

o 



.337 

o 



Ordering the elements in the following way {(0, 0, 0), (0, 1,0), 

(1, 0, 0), (1, 1, 0), (0, 0, 1), (0, 1, 1), (1, 0, 1), (1, 1, 1)}, the matrix Tx, = 





.113 




.544 
.337 
.448 
.011 



.456 


.439 




is an invariant part of the transition matrix Tx = 




.663 


.989 



2^11 T12 
Tx, 

Using the projection tt : — » Xi, -K{x,y,z) = {x,y,l); and the isomorphism 

p{x, y, 1) = (x, y), the network X is projected over the network Xi. Checking the 
arrows the projection 7f is a .5-homomorphism. 

5.1. Mathematical background. 

Theorem 5.4. If <j)i : Xi X2 is an ei-homomorphism, and (j)2 ■ X2 ^ X3 is 
another e2-homomorphism. Then (p — (j)2 o : Xi ^ X3 is an e-homomorphism. 
Therefore the Probabilistic Regulatory Networks with the homomorphisms of PRN 
form the category PRN. 

Proof. The Probabilistic Regulatory Networks with the PRN homomorphisms is a 
category if: the composition is an homomorphism, and satisfy the associativity law; 
and there exists an identity homomorphism for each PRN. 

(1) Let : Xx ^ X2 be an ei-homomorphism, and let (j)2 X2 ^ A3 be an 
e2-homomorphism. If qt, Qk and fj are functions in each PRN, and such that 
^1 o fj = g^o (pi and <f)2° gk = Qt° <p2, then we will prove that: <p o fj = qt o <P- In 
fact, 

{(j)2 O O fj =4>20 O fj) = 

(t>2 o {gk o <Pi) = (<A2 o gk) o(j)i = 
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{qt o 02) o 01 = gt o (02 o 0i)- 

(2)We want to prove that 

X(ife(0(u),0(w))) >xiciiu,v)). 

Suppose that x{ci{u, v)) = 1. Then, since 0i is an homomorphism of PRN, we have 
that 

X(dj(0i(u),0i(i;))) > x{ci{u,v)) 
which is 1. Since 02 is an homomorphism of PRN, we obtain that 

X(ife(0(«),0(w))) =x(ifc(02(0i(n)), 02(01 («)))) 
>x(cj(0i(w),(0i(t;)) = l. 

Therefore we obtain that 

X(i/c(02(0l(«)),02(0l(t^)))) =1. 

Then the composition of two PRN-homomorphisms is an homomorphism. 

(3) To verify the third condition for e-homomorphism, we do the foUowing. If 
p(0(wi), 0(u2)) > 1, with ui, U2 G ^1, then we need to prove that there exists an 
e such that 

|p(ui,«2) -p(0(wi),0(w2))| < e- 

In fact: 

\p{ui,U2) -_p(0(ui),0(m2))| = \p{ui,U2) - p(0i (ui ) , 0i (u2 ) ) + 
p(0l(wi), 01(^2)) -_p(02(0l(wi)), 02(01 (W2)))| 
< |p(ui,M2) -p(0l(wi),0l(U2))| + 
|p(0l(ui),0l(U2)) -p(02(0l(ui)),02(0l(W2)))| < Cl + £2 

\p{ui,U2) - p{(f>{ui), 0(u2))| < ei + e2 
because 0i and 02 are e-homomorphisms. 

The associativity and identity laws are easily checked, therefore our claim holds, 
and PRN is a category. □ 

It is clear that, the PRN with the homomorphism between them form a category 
that we will denote PTZM. The category PRN is a subcategory of PTZAf, since 
an homomorphism is not always an homomorphism for some e € M enough small. 
But, if we don't include the condition for e to be enough small, the two categories 
are the same, because always an homomorphism is an e-homomorphism for some 
e G R. 

Theorem 5.5. LetXixA:2 = {XixX2, H,E) be a product of PRN Xi = {Xi,F,C) 
and X2 = {X2,G, D). If 6i : X ^ Xi are two PRN-homomorphisms, then there 
exists an homomorphism S : X ^ Xi x X2, such that (j)ioS = Si for i = 1,2. That 
is, the following diagram commutes 

Xi XX2 

/ T \ 
Xi ^ X ^ X2 

This homomorphism is unique. 
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Proof. The function S : X ^ Xi x X2 is defined as follows S{x) = {Si{x), d2{x)), 
X G X. (5 is an homomorphism, in fact: 

(1) Let X = {X,L,P) be a PRN. Since 5i and S2 are homomorphism, for all 
function It & L there exist two functions fi £ F and gj G G, such that Siol^ = fio5\, 
and 62 °lt = gj ° ^2- Then for the function It there exists the function {fi,gj) that 
satisfies S o It = {fi, gj) o 5. 

{6olt){x) = 5{lt{x)) = {Si{lt{x)),52{lt{x))) = 

{MS,{x)),gj{S2{x))) = {{fi,gj)oS){x) 

(2) In order to prove x(eij((5(a;), (5(a;'))) > xiPhi^^ x')), suppose x{Ph{x, x')) = 1. 
Then lt{x) = x', and 5{x') = 6{lt{x)) = {fi, gj){S{x)) by part (1). Therefore 
x{eij{S{x) , S{x'))) = 1, and our claim holds. 

It is easy to check that 6 = Si, in fact 

0i((5(a;)) = (t)i{Si{x),52{x)) = Si{x), 
for all xgX. □ 

If ^i, i = 1, 2, are ej-homomorphism then 

max\p{x,x') - p{(l)i{S{x)), ^i{S{x')))\ < ei. 

But 

\pix,x') -piSix),Six'))+piSix),S{x'))-p{MSix)),MSix'm < 
\p{x,x')-p{S{x).S{x'))\ + 
\p{6ix),S{x'))-p{Mm),MS{x'm < ei. 

Therefore 

\p{x,x') -pi6ix),Six'))\ < ei - \pi6ix),6{x'))-p{Mm),MS{x'm 

\pix,x') - p{S{x),5{x'))\ < -Ii. 

Therefore 5 is an 6-homomorphism. So, the theorem holds for e-homomorphism. 

It is an imincdiato consequence the following resiilt, also is true for e-homomorpliisms. 

Theorem 5.6. Let Xi®X2 = {X1XX2, H,E) be a product of PRN Xi = {Xi,F,C) 
and X2 = {X2 ,G,D). If ji : Xi ^ X are two PRN-homomorphisms, then there 
exists an homomorphism 7 : Xi © X2 — » X , such that 7 o = 7^ for i = 1,2. That 
is, the following diagram commutes 

Xi®X2 

1-1 7 1'2 
/ i \ 

Xx^X^X2 

This homomorphism is unique. 

Theorem 5.7 (Fundamental Theorem). All reducible PRN is either a product of 
its non trivial sub-PRN or a subnetwork of this product. 

Proof. It is trivial by definition of Product and sub-PRN. □ 
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6. Conclusions 

The intersection, and the union of two sub-PRN is a sub-PRN, therefore the class 
of sub-PRN of a particular PRN is a lattice. Reduction mappings described in [7] 
and defined for PBN using the influence of a gene, for example Xn, on the predictor 
function /> to determine the selected predictor, can be extended to PRN. In order 
to extend this procedure to more than boolean functions, we use the polynomial 
description of genetic functions given in the partial derivative is the usual in 
calculus and all the concepts in |Zj can be using for PRN. Similarly our definition 
of projection, the reduction mappings are e-homomorphisms, and we can use for 
genes with more than two quantization, since this extension is not a trivial work 
we develop the theory and methods in 1 . 
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